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Abstract 

We present a simple theoretical framework to describe the thermal noise of a microscopic mechanical 
beam in a viscous fluid: we use the Sader approach to describe the effect of the surrounding fluid (added 
mass and viscous drag), and the fluctuation dissipation theorem for each flexural modes of the system to 
derive a general expression for the power spectrum density of fluctuations. This prediction is compared 
with an experimental measurement on a commercial atomic force microscopy cantilever in a frequency range 
covering the two first resonances. A very good agreement is found on the whole spectrum, with no adjustable 
parameters but the thickness of the cantilever. 

1 Introduction 

Cantilevers of micrometer size are nowadays present in many applications, ranging from chemical and biological 
sensors [1] to scanning probe microscopy [2]. In many cases, these tiny mechanical systems operate in a fluid 
environment (air or water for instance) which has a great influence on their dynamical behavior: the viscosity 
of the medium will broaden the structural resonances while the added mass due to the fluid moving along 
with the cantilever will shift their frequencies. A few theoretical models have been proposed to account for 
these effects [3-12], and validated experimentally [9-15]. Among the prediction of these approaches, the power 
spectrum of thermal noise induced fluctuations is of particular interest for its applications in atomic force 
microscopy (AFM) or micro-electromechanical systems (MEMS). We will focus here on the approach by Sader 
et al. [3,5,8], and derive a generic formula for the thermal noise using the fluctuation dissipation theorem for 
each mode, extending the work of Paul and Cross [4] in a simpler framework than Dorignac et al. [7]. This 
careful analysis is indeed incorrect in the original work of Sader et al., as noted by Paul and Cross [4], Dorignac 
et al. [7] and Sader et al. themselves (see ref. 44 of [5]). We will then compare the prediction of the model 
to a noise measurement on a commercial AFM cantilever in air: an excellent agreement is obtained for the 
full frequency behavior of the two first resonances, as well as at low frequencies and in the intermediate region 
between the two pics. 



2 Theory 

We will first recall the general theory to describe the motion of a cantilever immersed in a viscous fluid and 
subject to an arbitrary external forcing, using the same notations and approach as J. Sader [3], and then derive 
the expression for the thermal noise of the mechanical system. The cantilever is sketched on FiglT] Its length 
L is supposed to be much larger than its width 6, which itself is much larger than its thickness h. We will 
limit oursclf in this study to the flexural modes of the cantilever: the deformations are supposed to be only 
perpendicular to its length (along axis z of FiglT]) and uniform across its width. Theses deformation can thus 
be described by the deflection w{x,t), x being the spacial coordinate along the beam normalized by L, and t 
the time. The equation of describing the dynamics of w is [16]: 

Ei dMx,t) , dM^,t) 



Figure 1: Geometrical model for the cantilever: rectangular beam with length L, width b and thickness h. 
Deformation are supposed to be only along axis z and to depend only of time t and spacial coordinate x, they 
are described by the deflection w{x, t). 



where E is the Young modulus, / = bh^ /12 the second moment of inertia of the cantilever, the mass per unit 
length and F the external force per unit length. The boundary conditions for eq. [T] correspond to the classical 
clamped (x = 0) and free end {x ~ 1) conditions: 



w{x=0,t) = 
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In the Fourier space, we can rewrite eq. [l]as 

k d'^w{x\Lij) 
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dx^ 



mcui^'w{x\uj) = F{x\uj)L 
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with uj the pulsation in the Fourier space, rric = IJ-L the mass and k = 3EI/L^ the spring constant of the 
cantilever. Solving this equation in absence of external forces leads to the well known normal modes of oscillation 
of a clamped beam in vacuum with spatial profiles 



4>n{x) = (cosC„x — coshCnx) H — ; — " . " (siu CnX — smhCnx) 



where C„ is the n solution of equation 



sin Cn + sinh C„ 



1 + cos Cn cosh Cn = 



(5) 



(6) 



which leads to Ci = 1.875, C2 = 4.694, . . . , Cn ^ {n — l/2)n. The normals modes (pn{x) for n = 1, ... ,00 form 
an orthonormal basis of the fmictions of x in [0, 1] [16]: 



4>nix)(j)m{x)dx = Sn,: 



The pulsations of these resonant modes are given by the dispersion equations 



" k 
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When the cantilever is moving inside a fluid, the force acting on the cantilever can be decomposed into 2 
components; the force corresponding to the hydrodynamic load i^hydro(a;|w) due to the motion of the beam in 
the fluid, and the actual external driving Fdrivo(2^|w). Following Sader's approach [3], the hydrodynamic load 
can be approximated by 



^hydro(a;|w) = -pUj'^b'^T{uj)w{x\Lj) 



(9) 



where T{lu) is the hydrodynamic function corresponding to the rectangular cantilever beam and p the density 
of the fluid. The real part of F corresponds to the added mass due the fluid moving along with the cantilever 
during its motion (normalized to m/ = 7rpL5^/4, the mass of a cylinder of fluid of diameter b and length i), 
and the imaginary part F; to the viscous drag. To solve the equation of motion of the cantilever in the fluid, 
let us decompose the deflection w{x\uj) and external driving _Fdrivc(a;|^) on the orthonormal basis of the normal 
modes </>„: 

oo 

W{x\uj) = ^ Pn{uj)(f)n{x) (10) 

71=1 
^ OO 

{x\uj) = -^T]n{uj)(t)n{x) (11) 
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We note here that this decomposition is also valid in the real time space, with amplitude /3„(t) and external 
driving r]n{t) of each mode simply being the inverse Fourier transform of /3„(w) and rjn^uj). These 2 variables 
are coupled by the hamiltonian H of the system. Indeed, let us compute the infinitesimal work SW of i^drive 
when the deflection changes by 6w: 

dW = [ LdxFdrivc{x,t)Sw{x,t) (12) 



"'0 

oo 

= 5]77„(t)<5/3„(0 (13) 

ji=i 

For a reversible transformation, we can write dH = SW, hence 

'»» 

This last equation shows that the amplitude /3„ and external driving ?;„ of each mode are coupled variables 
by the hamiltonian of the system. We can thus apply the Fluctuation Dissipation Theorem [17]: the power 
spectrum density Sf3^ of fluctuation of amplitude for mode n is given by 

2kBT AjHx /-..x 
Sf}„M^ Im — ^ (15 

where fc^ is the Boltzmann constant, T the temperature of the thermostat and Im(.) is the imaginary part of 
its argument. 

We will now explicit the right part of eq. 1151 Using the orthonormalization of the base (/)„, we can easily 
show that eq. [4] leads to the following equation for each mode: 

|c,t - (m, + m/F(cj))w2^ /3„(a;) = 77„(a;) (16) 

The amplitude /3„ of each mode is thus governed by the equation of motion of an oscillator with stiffness 
mass rric + mfTr{uj) and damping coefficient m/wFi(a;), forced by external driving rjn- If F(a;) is a slow varying 
function of the frequency, close to its resonance each mode behaves like an harmonic oscillator. Using the above 
expression, eq. [15] can be rewritten as 

S. (w) ^ ^M^Ini(^ i — ) (17) 

'^"^ ' TTLO - (TO,. + TO^F(a;))t^2^ ^ ' 
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with T = mf/mc- In general, dissipation couples the modes calculated from the non-dissipative equations, but 
within the approximation of position independent dissipative forces (eqlH]), the noises are uncorrelated. We can 
thus write a generic expression of the power spectrum density of the deflection using eq. 1101 

oo 

sUx,uj) = Y,SpA^Mn{x)f (19) 

n=l 



This last equation is in agreement with the one derived by Dorignac et al. [7] in a more complex framework. 



3 Experiment 



To test the validity of the above formalism, we measure the thermal noise of a commercial AFM cantilever 
(Budget Sensor BS-Cont) with the following nominal geometry: L = 450 ^m, b = 50 /xm, h = 2 ^m. The 
cantilever is made of Silicium and is coating-less, thus for theoretical expressions we used the tabulated values 
of Si for the young Modulus {E = 169 GPa for Sino, the usual crystalline orientation along the length of AFM 
cantilever) and density {pc — 2340kg-m~^). The length and width were checked with an optical microscope, 
and the thickness is inferred from the resonance frequency of the cantilever (from eq. [S]), which eventually leads 
to L = 450 fim, 6 = 48 ^m, h = 2.07 fim. 

The measurement is performed with a home made intcrferomctric deflection sensor [18], inspired by the 
original design of Schonenbcrgcr [19] with a quadrature phase detection technique [20]: the interferences between 
the reference laser beam reflecting on the base of the cantilever and the sensing laser beam on the free end of 
the cantilever directly gives a calibrated measurement of the deflection w{x = with very high accuracy. 
We illustrate in Figl2]the performance of our detection with the power spectrum density (PSD) of a rigid mirror 
(in black): the background noise of our system is as low as 1.3 x 10~'^^m^/Hz, just 50% higher than the shot 
noise limit of our detection system. Typical noise of a measurement in a commercial AFM would be at least 
two orders of magnitude higher. 




Figure 2: Power spectrum density (PSD) of deflection. The PSD of thermal noise driven deflection at the free 
end of the cantilever (— ) is measured with an home made intcrferomctric setup in a IkHz to 100 kHz range. 
We clearly see the two first flexural modes of the mechanical beam, well above the background noise of the 
apparatus (— ) acquired on a rigid mirror. The difference between the measurement and the background noise 
spectrums (— ) gives a full view of the behavior of the thermal fluctuations of the cantilever in this frequency 
window. For comparison, the typical noise floor of a well tuned and calibrated commercial AFM is plotted (- -): 
only the resonances could be studied, the remaining part of the spectrum would be useless. 

We plot in Fig[2]the PSD of thermal noise driven deflection at the free end of the cantilever (in green). We 
clearly see in this spectrum the two flrst modes of oscillation of the mechanical beam, around /i = 14 kHz and 
/2 = 87 kHz, but we also have a lot of information out of these resonances as the measurement is always above 
the background noise of the system. If we subtract this background noise spectrum from the measurement 
(blue curve), we have an estimation of the thermal noise of the system on the whole 1 kHz to 100 kHz frequency 
interval, which we can compare to the theoretical expectation of the model presented in the previous section. 
The same measurement in a commercial AFM would yield much less information, as it would be limited to the 
resonances, and we would not have any information between them or at low frequency . 

In FiglSl we compare the measured PSD of fluctuations (in blue) with the prediction of 3 different models 
(in red) for the thermal noise of the cantilever, all based on Sader's approach [3]. The hydrodynamic function 
T{(jj) used for these predictions was taken from eq. 20, 21a and 21b of ref. [3], as it correspond to the flrst 
two resonances and the extension to higher order modes of ref. [8] is not necessary. The first model (FiglJl^a)) 
corresponds to eq. 29a of the original article of Sadcr [3], where the frequency dependence of the thermal 



noise forcing as been omitted. Both resonances are really well described with this model, but the ofF-resonance 
behavior is not: the low frequency noise measured is not flat as expected by the model, which also slightly 
underestimates the noise the intermediate region. 

The second model (FiglSjb)) corresponds to eq. 17a of ref. [8], a corrected expression of the thermal noise 
proposed by Green and Sader (see ref. 44 of [5] ) to account for frequency dependance of the dissipation following 
the approach of Paul and Cross [4]. Let us first note that this expression is dimensionally incorrect, as it is 
homogeneous to m^ and not to m^/Hz, so we had to rescale the formula by 2.28 x 10"'' Hz~^ to match the first 
resonance. However, this model perfectly describes the frequency dependance of the first mode, catching the 
slow variation of the noise in the low frequencies. The behavior of the second mode could also be well described, 
but it would need a higher rescaling factor. 

We came to the conclusion that the scaling factor is mode dependent, and conducted the analysis of the 
previous part to derive the correct expression. Indeed, it can easily be shown that the frequency behavior of 
each mode (eq. [T7| is exactly the same as in ref. [8], but the scaling factors arc different. We plot in (FiglSjc)) 
the prediction of eq. [17] and eq. I19|, and see that we match the whole power spectrum density of fluctuation of 
deflection: both resonances, low frequency and intermediate region are very well described by the model. We 
stress here that the only parameter we have adjusted in the model is the thickness of the cantilever /i, to match 
the actual frequencies of resonance. The Sader approach to describe the effect of the embedding fluid on the 
behavior of the cantilever turns out to be very good on the whole frequency range probed here. 

4 Conclusion 

In this paper, we have presented a simple theoretical framework to describe the thermal noise of a mechanical 
beam in a viscous fluid. Under the assumption of an infinitely thin and long cantilever, we used the Sader 
approach to describe the effect of the surrounding fluid (added mass and viscous drag). Using the fluctuation 
dissipation theorem for each flexural modes of the system, we derived a general expression for the power spectrum 
density of fluctuations as the sum of the contribution of the different modes. This prediction has been compared 
with an experimental measurement on a commercial AFM cantilever in a frequency range covering the 2 first 
resonances. A very good agreement has been found on the whole spectrum, with no adjustable parameters but 
the thickness of the cantilever. This analysis can be very useful in the AFM area, where the thermal motion of 
the cantilever is of great practical importance (as a lower bound to measurable forces or as a mechanical driving 
on its own) : it is not limited to a fit of resonances only, but gives an a priori knowledge of the cantilever noise 
on the whole frequency spectrum. 
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Figure 3: Power spectrum density (PSD) of thermal noise driven deflection at the free end of the cantilever: 
measurement ( — ) vs 3 different models ( — H — ). (a) prediction after [3], (b) prediction after [8], (c) prediction 
from eg . [T71 and eq. [iniof the present work. Insets present zooms around the resonances. 
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